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Distributed Parameters 


Having learned something about how we generate signals with bipolar and 
field effect transistors, we now turn our attention to the problem of getting 
those signals from one place to the next. Ever since Samuel Morse (and the 
founder of my alma mater, Ezra Cornell) demonstrated the first working 
telegraph, engineers and scientists have been working on the problem of 
describing and predicting how electrical signals behave as they travel down 
specific structures called transmission lines. 


Any electrical structure which carries a signal from one point to another can 
be considered a transmission line. Be it a long-haul coaxial cable used in 
the Internet, a twisted pair in a building as part of a local-area network, a 
cable connecting a PC to a printer, a bus layout on a motherboard, or a 
metallization layer on a integrated circuit, the fundamental behavior of all 
of these structures are described by the same basic equations. As computer 
switching speeds run into the 100s of MHz, into the GHz range, 
considerations of transmission line behavior are ever more critical, and 
become a more dominant force in the performance limitations of any 
system. 


For our initial purposes, we will introduce a "generic" transmission line 
[link], which will incorporate most (but not all) features of real transmission 
lines. We will then make some rather broad simplifications, which, while 
rendering our results less applicable to real-life situations, nevertheless 
greatly simplify the solutions, and lead us to insights that we can indeed 
apply to a broad range of situations. 

"Generic" Transmission Line 


The generic line consists of two conductors. We will suppose a potential 
difference V(x) exists between the two conductors, and that a current I(z) 
flows down one conductor, and returns via the other. For the time being, we 
will let the transmission line be "semi-infinite", which means we have 
access to the line at some point z, but the line then extends out in the x 
direction to infinity. (Such lines are a bit difficult to handle in the lab!) 


In order to be able to describe how V(a) and I(a) behave on this line, we 
have to make some kind of model of the electrical characteristics of the line 
itself. We can not just make up any model we want however; we have to 
base the model on physical realities. 


Let's start out by just considering one of the conductors and the physical 
effects of current flowing though that conductor. We know from freshman 
physics that a current flowing in a wire gives rise to a magnetic field, H 
({link]). Multiply H by p and you get B, the magnetic flux density, and 
then integrate B over a plane parallel to the wires and you get ®, the 
magnetic flux "linking" the circuit. This is shown in [link] for at least part 
of the surface. The definition of L, the inductance of a circuit element, is 
just 

Equation: 


where @ is the flux linking the circuit element, and I is the current flowing 
through it. Our only problem in finding ¢ is that the longer a section of wire 
we take, the more @ we have for the same J. Thus, we will introduce the 
concept of a distributed parameter. 


distributed parameter 
A distributed parameter is a parameter which is spread throughout a 
structure and is not confined to a lumped element such as a coil of 
wire. 


Example: 

For instance, we will hereby define E as the distributed inductance for 
the transmission line. It has units of Henrys/meter. If we have a length of 
transmission line zo meters long, and if that line has a distributed 
inductance of LZ H/m, then the inductance L of that length of line is just 
ie Lixo. 


Build Up of Magnetic Field 
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Likewise, if we have two conductors separated by some distance, and if 
there is a potential difference V between the conductors, then there must be 
some charge +(Q) on the two conductors which gives rise to that potential 
difference. We can imagine a linear charge distribution on the transmission 
line, p (C/m), where we have p Coulombs/m on one conductor, and —p 
Coulombs/m on the other conductor. For a line of length x9, we would have 
Q = +(pzo) on each section of wire. Whenever you have two charged 
conductors with a voltage difference between them, you can describe the 
ratio of the charge to the voltage as a capacitance. The two conductors 
would have a capacitance 

Equation: 
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and a distributed capacitance C’ (F/m) which is just ze A length of line x9 


long would have a capacitance C = C'xg Farads associated with it [link]. 
Find the Flux Linkage 
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Line Capacitance 


V(x) +Q 


Thus, we see that the transmission line has both a distributed inductance D 
and a distributed capacitance C’ which are tied up with each other. There is 
really no way in which we can separate one from the other. In other words, 
we can not have only the capacitance, or only the inductance, there will 
always be some of each associated with each section of line now matter 
how small or how big we make it. 


We are now ready to build our model. What we want to do is to come up 
with some arrangement of inductors and capacitors which will represent 
electrically, the properties of the distributed capacitance and inductance we 
discussed above. As a length of line gets longer, its capacitance increases, 
so we had better put the distributed capacitances in parallel with one 
another, since that is the way capacitors add up. Also, as the line gets 
longer, its total inductance increases, so we had better put the distributed 
inductances in series with one another, for that is the way inductances add 


up. [link] is a representation of the distributed inductance and capacitance 
of the generic transmission line. 
Distributed Parameter Model 


We break the line up into sections A(z) long, each one with an inductance 
EA(z) and a capacitance CA(z). If we halve A(x), we would halve the 
inductance and capacitance of each section, but we'd have twice as many of 
them per unit length. Duh! The point is no matter how fine we make 
CA(z), we still have Ls and Cs arranged like we see in [link], with the two 
kinds of components intermixed. 


We could make a more realistic model and realize that all real wires have 
series resistance associated with them and that whatever we use to keep the 
two conductors separated will have some leakage conductance associated it. 
To account for this we would introduce a series resistance R (ohms/unit 
length) and a series conductance G (ohms/unit length). One section of our 
line model then looks like [link]. 

Complete Distributed Model 


Although this is a more realistic model, it leads to much more complicated 
math. We will start out anyway, ignoring the series resistance R and the 
shunt conductance G.. This "approximation" turns out to be pretty good as 
long as either the line is not too long, or the frequencies of the signals we 


are sending down the line do not get too high. Without the series resistance 
or parallel conductance we have what is called an ideal lossless 
transmission line. 


Transmission Line Examples 


As an example, and also because it even has some practical importance, 
let's look at one kind of transmission line. It is called a stripline and it looks 
like [link]. It consists of a flat conductor, located between two ground 
planes. It is supported by an insulating dielectric with dielectric constant e. 
This is kind of like the situation you would find on a multi-level PC board, 
where perhaps the bus lines would be running on an inner layer with ground 
planes above and below them. 

A Stripline 
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Between the center conductor and the ground plane, there will be some 
capacitance, C’. If we can assume that the electric field is more or less 
confined to the regions between the strip conductor and the ground plane 
(which occurs when the ratio of i is not too small) then for either 
capacitor (assuming unit length into the picture) we will get a value 
Equation: 
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since the value of a capacitor is just the dielectric constant times the area of 
the plates, divided by the spacing of the plates. 


Looking quickly at [link] you might think the two capacitors are in series, 
but you would be wrong! Note that each capacitor has one lead connected 


to the center conductor and the other lead connected to ground, and so the 
two capacitors are in fact, in parallel, and hence their capacitances add. 
Thus, for the capacitance per unit length for this line, we can write: 
Equation: 
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It can be shown (although we won't do it here) that for any transmission 
line where the electric and magnetic fields are perpendicular to one another 
(called TEM or transverse electromagnetic) the speed of propagation of 
the wave down the line is just 

Equation: 
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Where é; is called the relative dielectric constant for the material. Well, 
we also know that 


Equation: 
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From which we can write 
Equation: 
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We can now insert this value for Z into the expression for Zp, the 
impedance of the line. 


Equation: 


And so, we have derived an equation for the impedance Zp of the line in 
terms of the dimensions W and B, the dielectric constant of the insulating 
material, ¢, and c, the speed of light. How good is this expression, and in 
particular how good is our assumption that the electric field is all confined 
to the region under the conductor? Not so great actually [link]. 

Exact and ae Impedance For a Stripline 
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[link] shows the results from using [link] and a more exact calculation, 
which takes into account the fringing fields. As you can see we have to get 
the ratio an up to about 4 or so before the two match. But at least we get 


the right behavior and we're not totally out of the ball park. 


Telegrapher's Equations 


Let's look at just one little section of the line, and define some voltages and 
currents [link]. 
Applying Kirchoff's Laws 


. ns dl (x,t) 
<a VL = Lh Ax at 
| (x,t) 
—s ——®_| (x +Ax,t) 
+ LAX bac Ay 2M Gee Ant) * 
ét 
V (x,t) C Ax 
V (x +Ax,t) 
—+—Ax 


For the section of line A() long, the voltage at its input is just V(a,¢) and 
the voltage at the output is V(a + A(a), t). Likewise, we have a current 
I(x, t) entering the section, and another current I(a + A(z), t) leaving the 
section of line. Note that both the voltage and the current are functions of 
time as well as position. 


The voltage drop across the inductor is just: 
Equation: 
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Likewise, the current flowing down through the capacitor is 
Equation: 
OV(x+ A(z), t 


Now we do a KVL around the outside of the section of line and we get 
Equation: 


V(a,t) —-V, —V(a2+ A(z), t) = 0 


Substituting [link] for Vz and taking it over to the RHS we have 
Equation: 


V(a,t) — V(x + A(a),t) = LA(z) 


Let's multiply by -1, and bring the A(z) over to the left hand side. 
Equation: 


V(a« + A(a),t) — V(z, t) OI (zx, t) 
A(z) --(4 ot ) 


We take the limit as A(x) — 0 and the LHS becomes a derivative: 
Equation: 


Now we can do a KCL at the node where the inductor and capacitor come 
together. 
Equation: 


OI(a + A(z),t) 


I(az,t) — CA(z) - 
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And upon rearrangement: 
Equation: 


I(x+A(x),t)-I(z,t) OV(« + A(z), t) 
A(z) 7 (c 


Now when we let A(x) — 0, the left hand side again becomes a derivative, 
and on the right hand side, V(a2 + A(x), t) — V(za, t), so we have: 
Equation: 


[link] and [link] are so important we will write them out again together: 
Equation: 
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Equation: 


These are called the telegrapher's equations and they are all we really 
need to derive how electrical signals behave as they move along on 
transmission lines. Note what they say. The first one says that at some point 
x along the line, the incremental voltage drop that we experience as we 
move down the line is just the distributed inductance L times the time 
derivative of the current flowing in the line at that point. The second 
equation simply tells us that the loss of current as we go down the line is 
proportional to the distributed capacitance C’ times the time rate of change 
of the voltage on the line. As you should be easily aware, what we have 
here are a pair of coupled linear differential equations in time and 
position for V(z,t) and I(z, t) 


Telegrapher's Equation in Real Lines 


So far, the transmission lines we have looked at have been "ideal". That is 
they have been lossless and dispersionless. Lest you leave the course with a 
false idea of how things really work, we should go back to our model and 
try to get things adjusted just a bit. 


As you can probably imagine, a real transmission line is going to have some 
series resistance, associated with the real losses in the copper wire. There 
may also be some shunt conductance, if the insulating material holding the 
two conductors has some leakage current. We will need to include these 
effects along with the distributed inductance and capacitance which we 
have already talked about. Fixing up the model accordingly, we now draw a 
section of line A(a) long as shown in [link]. Taking the same voltage loop 
and current sum that we did back in the discussion of transmission lines, we 
come up with the following version of the telegrapher's equations. 
Equation: 
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and 
Equation: 
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Real Line Diagram 
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Clearly, we would like to simplify things if we can. Let's again make a 
sinusoidal time excitation assumption, and let I(a,t) and V(a, t) become 
phasors. Since the time variation is now represented by a simple e’*”” the 
time derivatives become just 2w. We have 


Equation: 

ove) = —((R+iwL)I(z)) 
and 
Equation: 

on) = —((G + wC)V(2)) 


The way to get a solution is, of course, just like we have always done. Take 
the derivative with respect to x of [link] 


Equation: 
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and then plug in [link] 
Equation: 
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The obvious solution to this (See how easy this gets after you've done it 
once or twice) is 
Equation: 


V(x) = Voe*) 


with 
Equation: 


y= 4/(R+ iwL) (G+ iwCl) 


This number, + is called the complex propagation constant. Obviously, in 
general, it will have both a real and an imaginary part: 
Equation: 


y= a +48 


and we have 
Equation: 


V(a) = Voe ore 


Let's choose the minus sign in the exponent, and write the two terms as a 
product. 
Equation: 


We see we have something similar to what we had before, but with just a 
minor difference. The e~‘*4") term is the propagating term which tells us 
how the phase angle of the phasor changes as we move along the line, and 
acts just like the 6 term we had before. Thus 

Equation: 


and 
Equation: 
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The a is called the attenuation coefficient, and obviously, the e~ (@*) term 
in [link] causes the amplitude of the wave to decrease as it moves down the 
line. [link] is a sketch of what a wave would look like if it is both 
propagating down the transmission line and also being attenuated. In a 
distance the amplitude of the propagating wave has fallen to e~? of the 


value it had when it started. 
Wave Decay 
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Let's take the minus sign solution in [link] and substitute back into [link] 
Equation: 
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From which we get 
Equation: 
R+iwL)(G+iwC 
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Thus we can say 
Equation: 


V(a) = Zol (x) 


where 
Equation: 
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In general, in order to find a, 8, Ro , and Xo , we would have to find the 
square root given in [link] and [link] for specific values of R, L, G, and C’. 
On the other hand, we could maybe come up with some reasonable 
approximations which might suffice for cases of real interest. Obviously, if 
a line is very lossy, we would not be very interested in using it, and so 
except in some very special cases where an extremely lossy line is 
unavoidable (usually having to do with signals at very high frequencies) we 
might see if we can find a low loss approximation. 


A/C Line Behavior 


If we are going to try to use phasors on a transmission line, then we have to 
allow for spatial variation as well. This is simple to do, if we just let the 


phasor be a function of x, so we have V(x). How the phasor varies in x is 
one of the things we now have to find out. 


Let's start with the Telegrapher's Equations again. 


Equation: 
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For V(a, t) we can now substitute V(x)ei and for I(z, t) we plug in 


I(x)e“*. So we get: 


Equation: 
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We take the derivative with respect to time, which brings down a 2w and 
then we cancel the e“” from both sides of each equation: 


Equation: 


ove) —_ (iwLT(z)) 
and 
Equation: 

iG = — (iwCV(2)) 


Viola! In one simple motion, we have completely eliminated the time 
variable, t, from our equations! It is not really gone, of course, for once we 


figure out what V (z) is, we have to multiply it by e*“* and then take the 
real part before we can extract once again, the actual V(a, t) that we want. 
Nonetheless, insofar as the telegrapher's equations are concerned, t has 
disappeared from the radar screen. 


To solve these we do just as we did with the transient problem. We take a 


~ 


or) on the right 


hand side, for which we can substitute [link], which leaves us with 
Equation: 


derivative with respect to x of [link], which gives us a 
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(- times - is +, but 22 = —1 and so we have a - in front of the w?). We then 
re-write [link] as 
Equation: 
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The simplest solution to this equation is 
Equation: 


V(z) mn Voe~ (iwVLCz) 


from which we can then get the actual voltage signal 
Equation: 


V(a,t) = V(x)e™ 


= Vo a (wt+wVLCz) 


Note that we could factor out an eivV LC from the exponent, which, since it 
is just a constant, we could include in Vo (and call it Vj, switch the order of 
x and t, and write [link] as 

Equation: 


Viet) = Vieil®* zis!) 


which looks a lot like the "general" f(a + vt) solution we were talking 
about earlier! 


The number w+/ LC is special. It is usually represented with a Greek letter 
@ and is called the propagation coefficient. Thus we have 
Equation: 


V(a, t) _ Voetuttbe) 


As previously, a point on the wave of constant phase requires that the 
argument inside the parenthesis remains constant. Thus if V(a1, t1) is 
going to equal V(29, t2) (i.e. what was at point x at t; is now at £2 at time 
ty» it must be that 

Equation: 


wt, + Px, = wtp + Bae 


or 
Equation: 
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Which one again, defines the phase velocity of the wave. Other 
relationships to keep in mind are 
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The first comes from the fact that the wave varies in x as e***. Thus when 
x = ¥, the wavelength, Gy just increases by 27, to get the phasor to go 
through one full rotation. Note also, as before, the choice of the minus sign 
in the + in [link] represents a wave going in the x direction, while the 
choice of the + sign will give a wave going in the —z direction. Clearly, by 
starting out taking the x-derivative of the equation for I(x, t) we would end 
up with 

Equation: 


I(a, t) = ew) 


Let's consider the two phasors then, and define the voltage phasor 
associated with the positive going voltage wave as 
Equation: 
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Veunle) = V te U2) 


and the negative voltage phasor as 
Equation: 
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We should keep in mind that both V* and V~ can be, and probably are, 
complex numbers. (From now on we will drop the little ~ over the variables 
because its very tedious to get it to show up with this word processor. You 
will just have to keep in mind that any variable we do not explicitly put 
inside absolute value markers (i.e. |V *|) is going to be, in general, a 
complex number). We will, of course, have similar expressions for the 
positive and negative going current waves. 


Let's consider the positive going current and voltage waves, and plug them 
into [link]. 
Equation: 
+ e—(i8z) 
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The x-derivative brings down a — (i), the e~ (Pz)'s cancel, and we have 
Equation: 
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But, since 8 = w/ LC we have 


Equation: 
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as we had before. 


So, what has changed? Not much from the case of transients on a line. We 
will now assume we have a steady state problem. This means we turned on 
the generator a long time ago. We assume that it has been connected to the 
line long enough so that all transient behavior has died away, and that 
voltages and currents are not changing any more (except oscillating at 
frequency w, of course). 


If the line is semi-infinite (or matched with a load equal to Zo) [link] then it 
is pretty obvious that 
Equation: 
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where Z, is the source impedance, and V, is the source voltage phasor. 
A Wave On a Semi-Infinite Line 


Terminated Lines 


If all we did was launch signals down semi-infinite transmission lines, we 
would not get very much useful work done. We really need to have a finite 
length line, and put something at the end...like a termination. So let's take a 
look at a terminated transmission line [link]. 
At the Load End of the Line 
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The line has characteristic impedance Zp and we assume that it is 
terminated with a load resistor R,;,. If we have connected a source to the 
other end of the line, then we will have launched a voltage wave V," anda 


current wave [ - down the line. If the line is Z long, then it will take a time 


7 = =, where Up = —1_ for the signal to get to the end of the line. 
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What happens when the signal gets to the load? We can assume some 
voltage Vz, will appear across the load resistor, and hence a current Jy, will 
flow through it. The most logical thing to assume would be that V, = Ve. 
But, we quickly run into a contradiction. If Vy, = sae then 

Equation: 
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But Kirchoff says the sum of the currents into the output terminal must 
equal zero, thus [ ‘i must equal Jy. This then says 
Equation: 


y= & 


which can only be true if 2; = Zo, which while possible, will not be the 
case in general. 


What are we to do? We have an obvious contradiction. The telegrapher's 
equations permit two solutions to the transmission line problem: a signal 
going in the z direction, V,*, and a signal going in the —z direction, Vi 
The only way out of our problem is to assume that when the Aa signal gets 
to the load, a new signal, going in the —z direction, is created, which then 
heads back towards the load. So, let's putina V, and J, [link]. 

Reflected Waves are Generated 
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Now there is a little problem with signs here that we have to deal with. We 
can either draw J, so that it points in the —z direction, and say 
Equation: 


or we can draw J, going in the z direction, and say 
Equation: 


Since most of the time, it is better if we define all currents going in the 
same direction, and since this equation says that V~ and J” are related by 
— Z, the latter choice seems to be the better one. 


Going by that convention then, we equate voltage on either side of the 
termination 
Equation: 


and we sum currents into the output terminal 


Equation: 
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This makes 
Equation: 
Rr, = 
Ve (Vir -V, ) 


which we can substitute into [link] to get: 
Equation: 


or 


and this can be solved for V, as 
Equation: 


Rr — Zo 
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where I yz, is called the load voltage reflection coefficient. 


We could also have solved for J," in terms of J, and we would have found: 
Equation: 
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where J yy, is the load current reflection coefficient. Note that J, = —JIyz, 


Let's take a break from equation manipulating, and think about what we 
have here. First of all, although the result we have obtained is very 
important, the method we used to get there was even more so. What did we 
do? We postulated a voltage and current on the line, and then took a look to 
see if that solution resulted in a reasonable result. In this case it did not. We 
were in gross violation of Kirchoff Laws! We rescued ourselves by taking 
the only possible escape route: we added an additional voltage and current 


to the solution. Since V~ and J™ are related to each other in a different 
manner than V* and I* (by —Zp rather than Zo), this gave us an 
additional degree of freedom so that we could simultaneously satisfy both 
the transmission line I-V relationships as well as the load I-V relationship. 


Let's take a look at Z\;, for a minute. Over what range can it vary? A glance 
at [link] shows that it depends on the range of AR. If we exclude the 
possibility of negative resistance (a reasonable exclusion) then 

0< Rz, < ow. As R, varies over this range, the voltage reflection 
coefficient goes from —1 < Jy, < 1. When # is lessthan Zo, the 
reflection coefficient is negative. When Ry is greaterthan Zo, the 
reflection coefficient is positive. When Ry; = Zo the reflection coefficient 
is zero, and the line is said to be matched. In a matched transmission line, a 
signal traveling down it is completely absorbed by the load, and nothing is 
reflected from it. In this case we can have the incident voltage and current 
signals just equal load voltage and currents without the need to add 
reflected waves. For an unmatched transmission line, a signal incident on 
the load is (partially) reflected, and a new signal starts moving back down 
the line towards the source. An example of a transmission line is a buss ina 
computer. What would be the implications of unmatched terminations on 
various connections between the buss and the computer circuitry? 


Two special cases of terminated transmission lines that are of interest are 
Ry, = 0 (ashorted line) and Rz = oo (an open line). For Ry; = 0, 

I, = —1 and so a signal with the same amplitude but opposite sign as the 
incident wave, is reflected back down the line. For Ry; = oo, a signal with 
the same amplitude and sign gets reflected back down the line. It is easy to 
remember which is which, and also to make sure you have the right order in 
the equation for the reflection coefficient, when you keep in mind that the 
voltage across the load is the sum of V,* and V, . If the line is terminated 
with a short circuit, Vz must equal 0. Since 

Equation: 


VY, =I Vit 


and 
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yy = ——, 
Rr+Zo 
From these equations we can see that with R; = 0 (a short), I, does 
indeed —1, andsoV, = —-V," and hence V; = Ve +V, =Oasit 
should. 


Now, let's go back to our terminated transmission line. What is going to 
happen to V, after it leaves the load? It obviously travels back down the 
line towards the generator. What happens when it gets back there? Time to 
look at [link]. 
Back at the Generator 
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At the generator end of the line we can not leave out V,", as it is still here, 
so long as the source remains connected. Once it is launched, it stays on the 
line forever. We just have to add in the new V, . What happens now, if we 
test Kirchoff's voltage law here at the generator end? 

Equation: 


Verighg=V, HV S0 


Substituting [5 = I - + I, and then using the impedance relationship 
between the current and the voltages we get: 
Equation: 


F ve of vy = 
Note that again, we used ~~ for J," and —-y~ for I; . 


We know 
Equation: 


a 


= —f, 
: Z+Rs = 


which if we substitute into [link] and re-arrange a little bit we get 
Equation: 


Re Zo Re. 
ee ee 
s( FoR ae} Z } 1 


The stuff inside the big parentheses sums to zero (as in fact it should, this is 
just the solution to the initial Vi" generation problem) and we are left with 

the uncomfortable conclusion that V; must be zero! What are we going to 

do? We will just have to add another wave Vise to the solution! ({link]) 

Yet Another Wave is Formed 


Now we have two ways to proceed from here. The first would be the dumb 
way, and try to solve this whole problem. But then we could also be smart 
and note that all of the equations we have relating voltages and currents in 


this problem are linear equations and hence we can use the principle of 
superposition [link]. 
Using Superposition 

_!o—<—<—— 
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(a) (b) 


Thus, the top sketch in [link] can be broken into the two circuits (a) and (b) 
at the bottom of the figure. Then, everything is trivial. We have already 
solved (a) - that was just our initial Vi launching problem. Part (b) is just a 
reflection problem again, with a wave of V, incident on a load of value Rg 
through a transmission line with characteristic impedance Zo. The reflected 
wave in this instance is Vis which is why we had to use numbered 
subscripts in the first place. 


If you do not believe using superposition is valid, you can try doing the 
problem over again, but it should be pretty obvious that we can write 


Equation: 
Rs —Zo.,_ _ 
Me 9 ¥Y = Pav. 
2 Rg +Zo 3 me 


with 
Equation: 


Practical Filters 


An ideal filter simply removes all unwanted frequencies, preserving the 
remaining frequencies exactly. This would resemble some sort of a finite 
rectangle function in frequency. However, a simple, finite rectangle 
function in frequency is an infinite sinc function in time. This is a problem. 
A sinc function is an infinite length signal in both the positive and negative 
directions, making it impossible to create in the real-world. This leads us to 
as what would happen if we just made this sinc function causal by 
"chopping it off" somewhere. What we find when we do this is that the 
frequency domain representation is no longer a perfect rectangle: it now 
does not fall off immediately and shows some wiggling where it was flat 
before. 


A truncated sinc function ... leads to a smooth, wiggly, "rectangle" 
in time... in frequency. 
Ww 4) HU 
e, w 
The Bands 


An ideal filter has two types of bands: the stop band defines the region of 
frequencies that are eliminated by the filter, while the pass band defines the 
region of frequencies that the filter allows through. Practical filters add one 
more, the transition band. This is the area where the filter is moving 
between the stop band and pass band. 


Filter Design Specifications 


In our look at filter design specifications, we will use the example of a 
lowpass filter. The extension to the other kinds of filters should be fairly 
straightforward. [link] shows the parameters for a lowpass filter in the 
frequency domain. 


Hijed mse 
wigpl= 
instae woxes, 


From this image, the passband is the region from —w, to wp, the transition 
bands are the regions from —w, to —w, and from w, to w, while the stop 
band is the region less than —w, or greater than w,. 


In the figure above, e, and es represent the acceptable tolerance (or error) 
around the desired level that the passband and stopband respectively may 
vary. The behavior within the transition band is not specified, allowing 
anything to occur there, as long as the width is within specifications. 


Example: 
This example will look at a moving average system. 


Time domain The frequency domain representation of the 
representation of moving average system is a sinc function. 
; ine 
the moving H(iw) = ee) ene ) 


average. ws 


In L) | 4 Cjs)| 
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Some notes about this system: 


e It is lowpass 

e It has linear phase with jumps of 7 radians when the sinc function 
changes sign 

e The duration of the filter is inversely proportional to its bandwidth 

e This filter is finite impulse response (FIR) 

e It cannot be built with passive R, L, C circuits 

e We do not have independent control over all four design 
specifications 


We are going to design a moving average filter with the following design 
specs: WwW, = 1007, e, = 0.1, e, = 0.1 
lH (9 


With this specification, we are allowing w, to be a dependant variable 
(since we need one). We can now find the equation for this moving average 
system. 

We begin with 

Equation: 


aie 
CoN 12a) 
2 


Ww 


We will now solve for 7’ with 
Equation: 


= sin(50(z,T)) 


= 0:9 


For these specs, T’ ~ 0.007. This means that |H(iw)| does not stay below 
en— OO. Lintlo. = (iin. 
\HGr) | 


©.-\ 


O-\ 
yooTy TINT 


A graphical look at the transfer function of this lowpass 
filter with the passband and stopband noted. 


It is very clear from this representation that the transition band is huge ( 
6717). This is a very bad filter, especially when you consider that it cannot 
be implemented with passive circuitry. Fortunately better filters (e.g. 
Butterworth, Chebyshev and Elliptical) do exist. 


Beyond Lowpass Filter Design 


In the discussion of the different filters (Butterworth, Chebyshev and 
Elliptical) is common to see explanations based on lowpass filters. This 
explanation is very nice when first learning about them, because it is 
sufficient to understand the fundamentals of each of them. It is acceptable, 
because there exist fairly straightforward techniques to convert these 
lowpass filters into highpass, bandpass or bandstop filters. These techniques 


are the lowpass to highpass transformation, lowpass to bandpass 
transformation and lowpass to bandstop transformation. 


Butterworth Filters 


The Butterworth filter is a filter that can be constructed out of passive R, L, C circuits. The magnitude of the 
transfer function for this filter is 
Equation: 


Magnitude of Butterworth Filter Transfer Function 


where 7 is the order of the filter and w, is the cutoff frequency. The cutoff frequency is the frequency where the 


magnitude experiences a 3 dB dropoff (where |H(iw)| = A , 
T 
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frequency response |H(«)|? 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
frequency 


Three different orders of lowpass Butterworth analog filters: n = {1, 4, 10}. 
As n increases, the filter more closely approximates an ideal brickwall lowpass 
response. 


The important aspects of [link] are that it does not ripple in the passband or stopband as other filters tend to, and 
that the larger n, the sharper the cutoff (the smaller the transition band). 


This transfer function is often seen in its normalized form of 
Equation: 


Magnitude of Normalized Transfer Function for Lowpass Butterworth Filter 


1 


[H(is)| = —+ 
V14+ wn 


Butterworth filters give transfer functions (H (iw) and H(s)) that are rational functions. They also have only 
poles, resulting in a transfer function of the form 
Equation: 

1 


(s — 81) ($ — S2)--+ (8 — Sn) 


and a pole-zero plot of 


-1 -0.5 0 0.5 1 
Re(s) 


Poles of a 10th-order (n = 5 ) lowpass Butterworth filter. 


Note that the poles lie along a circle in the s-plane. 


Designing a Butterworth Filter 


Designing a Butterworth filter is a trivial task. Since we know that the filter contains only poles, we know that we 
can write it as 
Equation: 


1 
8” + ans" 14+ +--+ ajs+1 


H(s) = 


From this, we may look up the a; from a table (like the one below) for any desired n. We can also find them in 
Matlab by using the but tap command. The real challenge of designing a Butterworth filter comes with figuring 
out the optimal characteristics for the given application. 


n ay a2 a3 a4 a5 a6 a7 
2 1.414214 


3 2.000000 2.000000 


n at 


4 2.613126 
5 3.236068 
6 3.863703 
7 4.493959 
8 5.125831 
9 5.758770 


10 6.392453 


Exercise: 


Problem: 


a2 


3.414214 


5.236068 


7.464102 


10.097835 


13.137071 


16.581719 


20.431729 


a3 


2.613126 


5.236068 


9.141620 


14.591794 


21.846151 


31.163437 


42.802061 


a4 


3.236068 


7.464102 


14.591794 


25.688356 


41.986386 


64.882396 


a5 


3.863703 


10.097835 


21.846151 


41.986386 


74233429 


a6 


4.493959 


13.137071 


31.163437 


64.882396 


az 


5.125831 


16.581719 


42.802061 


Design a Butterworth filter with a passband gain between 1 and 0.891 (-1 dB gain) for0 < w < 10 anda 
stopband not to exceed 0.0316 (-30 dB gain) for w > 20. 


Solution: 


The first step is to determine n. To do this, we must solve for n using the passband and stopband criteria. We 
begin by finding the equation for the gain in the passband in dB, 


Equation: 


and for the stopband in dB, 


Equation: 


— 


Gp = 20log|H(iw)| 


2s 108g (1 2 (#)”) 


G, = 20log|H(iw)| 


== = yolieg (1 + (2)”) 


these equations can also take the form 


Equation: 


none el & 
(=) =10% -1 
We 


In this form, we may divide the passband equation by the stopband equation to get rid of the w,. From there, 
we can solve for n to get 


Equation: 


By plugging in, we find n = 5.9569. However, since n must be an integer, we round this up ton = 6 


The next step is to find w,. We can do this by substituting nm = 6 into the equations for the passband and 
stopband and solving for w,. This yields w. = 11.1919 for the passband equation and w, = 11.2478 for the 
stopband equation. The difference in these solutions is a result of n needing to be an integer. If we choose the 
solution from the passband equation, the passband will meet its requirements exactly, and the stopband will 
surpass its requirements. If we choose the solution from the stopband equation instead, the stopband 
requirements will be met exactly, while we will exceed the passband requirements. Therefore, we may choose 
either value or any value in between. For this example, we will choose w, = 11.2478. 


Now, we can find the normalized transfer function. Since we know this to be a sixth-order Butterworth, we 
can determine from the table that 
Equation: 
1 
s® + 3.863703s° + 7.464102s4 + 9.141620s3 + 7.464102s2 + 3.863703s +1 


H(s) = 


Finally, we can determine the final transfer function. 
Equation: 


1 


H(s) = 


° + 3.863703(=-2..)° + 7.464102(=2..)* + 9.141620(=3..)* + 7.464102( 


8 8 : 
( 11.2478 ) 11.2478 11.2478 11.2478 11.2478 ) 


Rather than multiplying this out and factoring, we will leave it in this form for readability, since the numbers 
can get quite large otherwise. 


Chebyshev Filters 


This module will cover Chebyshev filters with an assumed understanding of 
Butterworth filters. This module will also examine the lowpass example of 
these filters, leaving conversion to other types of filters for the Lowpass to 
Highpass Transformation, Lowpass to Bandpass ‘Transformation and 
Lowpass to Bandstop Transformation modules. 


Like Butterworth filters, Chebyshev filters contain only poles. However, 
while the poles of the Butterworth filter lie on a circle in the s-plane, those 
of the Chebyshev filter lie on an ellipse. 

jw 


A pole-zero plot of a lowpass 
Chebyshev filter The poles are 
equally spaced around an ellipse in 
the left half of the complex plane. 


The result of this repositioning of poles is a "rippling effect" in the 
passband of the magnitude of the transfer function. Since each local 
maximum in this rippling reaches the same value and each local minimum 
reaches the same value, this rippling is described as equal ripple. It is 
important to notice that there is no rippling in the stopband and to be aware 
that this transition band will be narrower than a comparable Butterworth 
filter. 


3 
? 
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A sketch of the magnitude of a lowpass 

Chebyshev filter. Notice that the "equal 

rippling" in the passband and not in the 
stopband. 


The magnitude of the transfer function of a Chebyshev filter takes the form 
Equation: 


| H(iw)| = a 
V1 + €2C;,?(w) 


where the polynomial, C,,(w) = cos(n arccos(w)) is known as the 
Chebyshev polynomial. 


The design of Chebyshev filters is generally done the same way 
Butterworth filters are (with tables and Matlab). The most relevant Matlab 
command is cheb1ap, but there are others. 


